(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



CORRECTED VERSION 



(19) World Intellectual Property Organization 

International Bureau 

(43) International Publication Date 
19 September 2002 (19,09.2002) 



mm 

MM 



PCT 




mil 




(10) International Publication Number 

WO 02/073239 Al 



(51) International Patent Classification 7 : G01V 1/16, 1/38 

(21) International Application Number: PCT/IB02/0L897 

(22) International Filing Date: 12 March 2002 (12.03.2002) 



(25) Filing Language: 

(26) Publication Language: 



English 
English 



(30) Priority Data: 

0106091.2 



13 March 2001 (13,03.2001) GB 



(71) Applicant (for all designated States except CA, FR, US): 
WESTERNGECO SEISMIC HOLDINGS LIMITED 

I — / — ]; P.O. Box 71, Craigmuire Chambers, Road Town, 
Tortola, British Virgin Islands (VG). 

(71) Applicant (for CA only): WESTERNGECO CANADA 
LIMITED [CA/CA]; 24th Floor, Monenco Place, 801 6th 
Avenue, S W, Calgary, Alberta T2P3W2 (CA). 



(71) Applicant (for FR only): SERVICES PETROLIERS 
SCHLUMBERGER [FR/FR]; 42, rue Saint-Dominique, 
F-75007 Paris (FR). 

(72) Inventors; and 

(75) Inventors/Applicants (for US only): ROBERTSSON, 

Johan [SE/NO]; Slasjonsveien 28A, N-0772 Oslo (NO). 
MUUS, Remco {NL/CH]; institute of Geophysics, 
ETH-Hoenggerberg, CH-8093 Zurich (CH). 

(74) Agent: SUCKLING, Andrew, Michael; Marks & Clerk, 
4220 Nash Court, Oxford Business Park South, Oxford 
OX4 2RU (GB). 

(81) Designated States (national): AE, AG, AL, AM, AT, AU, 
AZ, BA, BB, BG, BR, BY, BZ, CA, CH T CN, CO, CR, CU, 
CZ, DE, DK, DM, DZ, EC, EE, ES, FI, GB, GD, GE, GH, 
GM, HR, HU, ID, IL, IN, IS, JP, KE, KG, KP, KR, KZ, LC, 
LK, LR, LS, LT, LU, LV, MA, MD, MG, MK, MN, MW, 
MX, MZ. NO, NZ. OM, PH, PL, PT, RO, RU, SD, SE, SG, 
SI, SK, SL, TJ, TM, TN, TR, TT, TZ, UA, UG, US, UZ, 
VN, YU, ZA, ZM, ZW. 

[Continued on next page] 



(54) Title: METHOD AND APPARATUS FOR DETERMINING THE ORIENTATION OF A SEISMIC RECEIVER DURING A 
SEISMIC SURVEY 




13 



< 
OS 

m 




O 



(57) Abstract: The invention provides a method of determining the orientation of a seismic receiver (4,6) from seismic data acquired 
at the receiver (4,6). The determined orientation can be Laken into account in subsequent analysis of the seismic data. This avoids 
inaccuracies that can occur if the receiver orientation is estimated. In a preferred embodiment the horizontal spatial derivatives of 
the pressure measured at the receiver are used in the determination of the orientation of the receiver These may be calculated on 
either the source side or the receiver side. These arc combined with horizontal components of the particle displacement, velocity 
or acceleration (or a higher time-derivative of the particle displacement). Alternatively, horizontal spatial derivatives of the particle 
displacement measured at the receiver may be used in place of the horizontal spatiaJ derivatives of the pressure. The invention also 
provides a seismic receiver (6) having three or more closely-spaced, non-collinear pressure sensors (8,9,10). The provision of three 
or more pressure sensors enables the horizontal spatial derivatives of the pressure to be determined accurately. 
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(57) Abstract: The invention provides a method of determining the orientation of a seismic receiver (4,6) from seismic data acquired 
at the receiver (4,6). The determined orientation can be taken into account in subsequent analysis of the seismic data. This avoids 
inaccuracies that can occur if the receiver orientation is estimated. In a preferred embodiment the horizontal spatial derivatives of 
the pressure measured at the receiver are used in the determination of the orientation of the receiver. These may be calculated on 
cither the source side or the receiver side. These are combined with horizontal components of the particle displacement, velocity 
or acceleration (or a higher time-derivative of the particle displacement). Alternatively, horizontal spatial derivatives of the particle 
displacement measured at the receiver may be used in place of the horizontal spatial derivatives of the pressure- The invention also 
provides a seismic receiver (6) having three or more closely-spaced, non-collinear pressure sensors (8,9,10). The provision of three 
or more pressure sensors enables the horizontal spatial derivatives of the pressure to be determined accurately. 
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METHOD AND APPARATUS FOR DETERMINING THE ORIENTATION OF A SEISMIC RECEIVER 
DURING A SEISMIC SURVEY 

The present invention relates to a method of detennining the orientation of a seismic 
receiver, in particular to a method of determining the heading (that is, the orientation in 
a horizontal plane) of a multi-component seismic receiver. The invention can be 
applied to, for example, a receiver disposed on the earth's surface, to a receiver being 
towed through water, or to a receiver disposed within a borehole. The term "earth's 
surface" as used herein includes the sea-bed, land, and the transition zone. The 
invention also relates to a seismic receiver, and to a method of seismic surveying using 
such a seismic receiver. 

A seismic receiver generally comprises one or more seismic sensing elements disposed 
within a housing. When the receiver is deployed in a desired survey location by being 
disposed on the earth's surface, the coupling of the seismic sensing element(s) to the 
earth is provided by the housing of the receiver; the housing also provides physical 
protection for the sensing element(s). Receivers of this general type are generally used 
by attaching the receiver housings to a support cable at intervals along the length of the 
cable. The support cable is provided with electrical leads to enable output signals from 
the receivers to be transmitted to monitoring and/or recording equipment, and to enable 
power to be provided to the receivers. For a land-based survey, the cable is disposed on 
the land so that the receivers are located at their desired positions, In the case of sea- 
bed seismic data acquisition, the cable is lowered onto the sea-bed to deploy the 
receivers at their desired locations on the sea-bed. Alternatively, in the case of a towed 
marine receiver array, the cable is suspended at a desired depth below the water surface 
and is towed through the water, for example by a survey vessel. As a further 
alternative, the cable may be deployed within a borehole to provide a VSP (vertical 
seismic profile) receiver array. Receivers may also be deployed without being attached 
to a cable. For instance, autonomous sensors can record and store data locally, for 
example on a hard disk, for later retrieval, or can transmit data, for instance via radio, to 
another location for recordal there. 
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A seismic receiver provided with two or more seismic sensing elements is able to 
measure two or more parameters of the received seismic energy, and is thus known as a 
multi-component seismic receiver. One application of multi-component seismic 
receivers is as sea-bed seismic receivers - that is, as receivers intended to be disposed on 
the sea-bed. Sea-bed seismic receivers generally record the pressure and the elastic 
wavefield of the seismic energy incident on the receiver. The pressure is a scalar 
quantity, whereas the elastic wavefield is a vector quantity and it is therefore necessary 
to measure the components of the elastic wavefield in three non-coplanar directions. 
This is done using a three-component geophone that can measure the components of the 
wavefield in three non-coplanar directions. The three directions chosen are usually the 
x-direction, the y-direction, and the z-direction (vertical). In total, therefore, four 
components of the seismic data are measured. Four-component seismic data recording 
at the sea-bed has proven to be a very successful method for imaging through gas 

» 

saturated overburdens and for characterising hydrocarbon reservoirs through lithology 
identification and fluid discrimination. The multi-component data describing the elastic 
wavefield are especially useful, since they enable the separation of the compressional P- 
waves from the shear S-waves. 

In a conventional seabed acquisition system the receivers are firmly coupled to the 
cable, so that the orientation of each receiver, relative to the cable, is known. Usually, 
each receiver is aligned on the cable so that one sensing element records the component 
of the elastic wavefield along the cable (this is also known as the in-line direction, and 
is generally taken to be the x-direction), and another records the component of the 
wavefield transverse to the cable (this is also known as the cross-line direction, and is 
generally taken to be the y-direction). The third sensing element records the vertical 
component of the wavefield. 

In a conventional receiver array, the orientation of each receiver on the cable is defined 
when the receiver is attached to the cable, so that the receiver heading is defined before 
the cable is deployed. Thus, the heading of each receiver can in principle be determined 
from its position on the cable. However, twisting and coiling of the support cable may 
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occur as the cable is deployed, and this would mean that actual heading of a receiver 
may be different from the theoretical heading determined from the receiver's position 
on the cable. Furthermore, new seabed acquisition systems are being designed in which 
the receivers are decoupled from the cable. In these new acquisition systems, although 
each receiver will still measure the components of the wavefield along its internal x- 
and y-axes, knowledge of the receiver's position on the cable will not be sufficient to 
determine the orientation of the receiver's x- and y- axes. Thus, determining the 
receiver's heading from its position on the cable may well be inaccurate, with a possibly 
severe effect on the quality of the results of processing the data. 

Determining the receiver heading solely from the receiver's position is not possible in 
the case of an autonomous receiver. 

A first aspect of the present invention provides a method of determining the orientation 
of a seismic receiver, the method comprising: acquiring seismic or acoustic data at a 
receiver, the data including at least one of pressure and the particle displacement or a 
time-derivative of the particle displacement; and determining the orientation of the 
receiver from at least one of the pressure measured at the receiver and the particle 
displacement or time-derivative of the particle displacement measured at the receiver. . 

The present invention thus eliminates the need to estimate the heading of a receiver 
from its position on a support cable. Instead, the heading can be determined from 
seismic or acoustic data acquired by the receiver, and this allows more accurate 
determination of the heading of a seismic receiver. The invention enables the receiver 
heading to be determined even in the absence of any prior knowledge about the heading 
of the receiver, and is of particular benefit when applied to an acquisition system in 
which the receivers are decoupled from a support cable or in which autonomous 
- receivers are used. 

Once the heading of the receiver has been determined, seismic data acquired by the 
receiver can be processed in any conventional way. The processing can take account of 
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the determined heading of the receiver, thereby improving the accuracy of the data 
processing. 

In the method of the invention, the heading of the receiver is determined from at least 
one of the pressure measured at the receiver and the particle displacement, or a time- 
derivative thereof measured at the receiver. The most convenient time-derivatives of 
the particle displacement to use are the particle velocity (i.e. the first time-derivative) 
and the particle acceleration (i.e. the second time-derivative) but, in principle, higher 
order time-derivatives could be used. 



In a preferred embodiment the spatial derivatives, in first and second different 
horizontal directions, of the pressure measured at the receiver are used in the 
determination of the receiver orientation. Components of the particle displacement, or 
particle velocity^ or particle acceleration, measured at the receiver in the first and 

■ 

second horizontal directions may also be used in the determination of the receiver 
orientation. The spatial derivatives, in the first and second horizontal directions, of the 
particle displacement measured at the receiver may be used instead of the spatial 
derivatives of the pressure. 



Faithful recording of the pressure and elastic wavefield is essential for a full and correct 
analysis of the seismic data obtained in a multi-component seabed seismic survey. 
Acquired seismic data can receive the full benefit of vector processing techniques, such 
as wavefield decomposition schemes, only if they provide a good vector representation 
of the wavefield. In practice, variations in geophone and hydrophone coupling, 
resonances in the receiver cable, and incorrect receiver orientation can affect the 
recordings. 



One particular problem that occurs in processing data acquired by a multi-component 
seismic sensor is that of vector infidelity. By "vector infidelity" is meant that one or 
more components of the recorded wavefield is/are distorted compared to the true 
particle motion. Sub-optimal design of the sensor housing or of the cable can cause 
such distortions, which tend to be particularly severe for specific components of the 
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recorded data. In particular, infidelity between the in-line (x-direction) and cross-line 
(y-direction) components of the elastic wavefield can be a serious problem. 

Various techniques for correcting for vector infidelity have been proposed. For 
example, J.E. Gaiser, in "Compensating OBC data for variations in geophone 
coupling", 68* Ann. Internal. Mtg, Soc. Expl. Geophys. ppl429-1432 (1998), and C. 
Bagaini et al, in "Assessment and calibration of horizontal geophone fidelity in seabed- 
4C using shear waves", 62nd EAGE Conference Glasgow, Extended Abstracts, paper 
L02 (2000), have proposed techniques for calibrating the recordings of different 
components of the wavefield in order to correct for vector infidelity. These techniques 
rely upon the minimisation of transverse energy, assuming all energy to travel in the 
radial-vertical plane. Both schemes, however, define the radial direction as the 
direction from source to receiver, which implicitly assumes that the earth is laterally 
invariant This assumption means that any subsurface structures present will cause 
systematic uncertainties in the calibration filters, as pointed out by Gaiser (supra). 

♦ 

The errors in the calibration filters have hitherto been minimised statistically, but tins 
can be done only if the acquired data has a good azimuth coverage. The present 
invention, however, makes it possible to perform a correct receiver calibration without 
the assumption of a 1-D earth model, for any number of shots with an arbitrary azimuth 
. coverage. 

A second aspect of the present invention provides a seismic receiver comprising three or 
more pressure sensing elements, the positions of the pressure sensing elements not 
being collinear. Such a receiver can provide information about the spatial derivatives of 
the pressure, and this information can be used in a determination of the heading of the 
receiver. Seismic data acquired by a receiver according to the second aspect of the 
invention can be processed by a method according to the first aspect of the invention to 
determine the receiver orientation. 

A third aspect of the present invention provides a method of seismic surveying 
comprising: emitting seismic energy at a first location; and acquiring seismic data at a 
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seismic receiver of the type defined above disposed at a second location different from 
the first location. 

Some or all of the seismic data acquired by a seismic surveying method of the third 
aspect of the invention can be processed according to a method of the first aspect of the 
invention so as to determine the orientation, in particular the heading, of the receiver 
used to acquire the data. Subsequent processing of the acquired seismic data can take 
account of the determined receiver heading. 

Other preferred features of the present invention are set out in the dependent claims. 

Preferred embodiments of the invention will now be described by way of illustrative 
example with reference to the accompanying Figures in which: 

■ 

Figure 1 is a schematic plan view of a first marine seismic surveying arrangement 
suitable for carrying out a method of the invention; 

Figure 2 is a schematic plan view of a second marine seismic surveying arrangement 
suitable for carrying out a method of the invention; 

Figure 3(a) illustrates a first energy function obtained for noise-free synthetic seismic 
data; 

■ 

Figure 3(b) illustrates a second energy function for the data of Figure 3(a); 

Figure 3(c) corresponds to Figure 3(a) but is for synthetic seismic data having a signal- 
to-noise ratio of 20dB; 

Figure 3(d) illustrates the second energy function for the data of Figure 3(c); 

Figure 3(e) corresponds to Figure 3(a) but is for synthetic seismic data having a signal- 
to-noise ratio of lOdB; 



WO 02/073239 



PCT/IB02/01897 



7 

Figure 3(f) illustrates the second energy function for the data of Figure 3(e); and 

Figure 4 is a schematic illustration of a multi-component seismic receiver according to 
an embodiment of the present invention. 

In the method of the present invention the receiver orientation in a horizontal plane (this 
is known as the "receiver heading**)* is determined from seismic or acoustic data 
acquired by the receiver. Information about at least one of the pressure at the receiver 
and the particle displacement, or a time derivative thereof, at the receiver is obtained 
from the data, and the receiver heading is determined from this information. In a 
preferred embodiment of the invention the spatial derivatives in two different horizontal 
directions of the pressure at the receiver are used in the determination of the receiver 
heading. The two directions are preferably perpendicular, and it will be assumed that 
the two directions are perpendicular in the following description of an embodiment of 
the invention. 

Figure 1 is a schematic plan view of a marine seismic survey. A seismic source 1 is 
being towed by a survey vessel 2. As the seismic source is towed through the water it is 
actuated to emit a pulse of seismic energy at each shot point 3. The distance between 
successive shot points 3 is equal to the product of the speed of the vessel 2 and the time 
interval between successive actuations of the source 1 . The direction of travel of the 
vessel 2 is defined to be the x-direction. 

* 

The seismic energy emitted by the seismic source 1 is received at a multi-component 
seismic receiver 4 mounted on a support cable 5. The receiver 4 contains sensing 
elements for sensing the components of the wavefield in two orthogonal horizontal 
directions. These directions are the receiver's internal x- and y-axes, and will be 
referred to as the directions x rcc and y^c respectively. The receiver 4 also contains a 
pressure sensing element. 
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The orientation of the receiver's internal x- and y-axes x Tec , y Tec relative to the x- and y- 
axes defined by the direction of movement of the boat (hereinafter referred to as the 
"global x-axis" and "global y-axis") is unknown. The angle between the receiver x-axis 
Xrtc and the global x-axis is denoted by a. 



The direction of propagation, or azimuth, of an incoming wave in the horizontal plane 
can be obtained from the horizontal spatial derivatives of either pressure P or 
displacement u. The angle between the azimuth of the wave and the positive global x- 
direction is shown as 9 in Figure 1. The azimuth can be determined from: 



ton^^f^* ' (1) 
d y P d y u t 



In equation (1), the notation ft denotes differentiation with respect to the variable *i\ 
and Ui denotes the i th component of the displacement (i = x, y or z). 

As may be seen from Figure 1, the azimuth of the wave can be determined, in principle, 
from spatial derivatives of either the pressure or the displacement measured at the 
receiver 4. However, the pressure recorded at a receiver is independent of the receiver 

r 

heading, since the pressure is a scalar quantity. Moreover, the measured pressure is 
generally less susceptible to variations from receiver-to-receiver in the coupling 
between a receiver and the elastic medium (e.g. the sea-floor or the ground) below or 
around the receiver than is the measured particle displacement. In preferred 
embodiments of the invention, therefore, the orientation of the receiver is determined 
from derivatives of the pressure rather than from derivatives of the displacement. 

* 

In the absence of seismic energy polarised out of the horizontal plane (this may arise 
owing to Sh-energy or shear wave splitting), it is the case that 

d x u y - d y u x = 0, (2) ~ 



and using the following expressions for the components of the horizontal slowness p 
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in the x- and y-directions, respectively: 



it is possible relate the horizontal spatial derivatives of pressure to the horizontal 
components of particle velocity, as follows: 



i-M. (4) 



Equation (4) allows the sensor heading to be determined. It can be re-written as: 



v x dyP-v y d x P = 0 (5) 



The receiver 4 will record the components of the particle velocity along the x^c- and 
yrec-axes, and will also record the pressure. Moreover the seismic surveying 
arrangement allows the spatial derivatives d^cP and d ync P of the pressure to be 
determined, as will be discussed below. If the quantity v mc dyrecP — v yre c dxrecP is 
computed, it should be zero if the receiver x- and y-axes are aligned with the global x- 
and y-axes but will be non-zero if the receiver x- and y-axes are not exactly aligned with 
the global x- and y-axes. 



To determine the heading of the receiver, therefore, the quantity v xrec dy re cP - ^yrec dxrecP 
is computed. If the quantity dy rec P - v yrec dxrecP should be zero this indicates that 
the receiver x- and y-axes are aligned with the global x- and y-axes, and that the angle a 
is equal to zero. If the initial computation of the quantity v wc d yre cP - Vync dxrecP should 
be non-zero, however, this shows that the receiver axes are not in alignment with the 
global x- and y- axes. In this case, dyrecP - dxrecP is recomputed for a new set 
of receiver x^ c - and y rec -axes for which the angle a has some non-zero value, and this 
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process is repeated until a value of the angle a is found that gives a value for the 
quantity v xrec d ynsc P — Vyrec dxreeP of zero, or acceptably close to zero. 



Since it is not known whether the initial calculation of - Vyrec dxre<P will give 

a positive or negative value, a preferred way of determining the receiver heading is to 
find the value of the angle a that minimises the cost function, or energy function, 
E\ = 6w dyncP - Vymc dxncFf. E\ will never be negative, and this simplifies iterative 
techniques for determining the angle a. 



Once the value of the angle a has been determined, the spatial derivatives of the 
pressure along the global x- and y-axes can be computed from the pressure derivatives 
along the receiver's x- and y-axes. The azimuth 8 of the wave can then be determined 
using equation (1). 

The above method can be carried out using any suitable computational technique. 

The above description of the method has assumed that the receiver is oriented so that its 
z-axis is vertical. In many cases this assumption will be correct, because deployment of 
the cable is more likely to cause errors in the orientation of the receiver's x- and y-axes 
than in the orientation of the receiver's z-axis. Furthermore, many seismic receivers 
now incorporate sensing elements mounted on gimbals so that correct orientation of the 
z-axis of the receiver is ensured. If, however, it is known or suspected that the 
receiver's z-axis may not be vertical it is possible to pre-process the data acquired by 
the receiver to correct for the mis-orientation of the z-axis before applying the method 
of the invention to determine the sensor heading. For example, any projection of a 
vertical component onto the x^ c and y^c components of data recorded by the receiver 
can be removed by calibrating the recorded data against the equation of motion using 
the technique disclosed by K. M. Schalkwijk, C. P. A. Wapenaar, and D. J. Verschuur, 
in "Application of two-step decomposition to multi-component ocean-bottom data: 
Theory and case study," J. Seism, ExpL, Vol 8, pp261-278 (1999). 
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As noted above, the derivation of equation (4) assumed that no seismic energy was 
polarised out of the horizontal plane. The determination of the receiver heading is 
therefore preferably performed on a portion of the seismic data for which this 
assumption is reasonable, by selecting a time-window of the seismic data in which the 
received seismic energy is unlikely to contain seismic energy polarised out of the 
horizontal plane. This can be done, for example, by selecting a time window which 
contains the direct wave (the wave that travels direct from the source 1 to the receiver 
4). If it is assumed that the sea floor in the vicinity of the survey area is flat and that 
there is no azimuthal anisotropy, selecting a time window that contains water layer 
multiple reflections (that is, events that involve reflection at the sea surface and/or the 
sea-bed) should also ensure that the data used in the determination of the receiver 
heading is unlikely to contain seismic energy polarised out of the horizontal plane. 

Once the orientation of the receiver has been determined by processing a selected 
portion of the seismic data, conventional processing steps can be applied to the seismic 
data. These subsequent processing steps are able to take account of the determined 
orientation of the receiver. Subsequent processing steps can be applied to the entire 
seismic data, or to one or more selected portions of the seismic data. 

In the above description, the receiver heading has been determined from seismic data 
acquired at the receiver. The invention is not limited to detennining the receiver 
heading from seismic data, and the receiver heading may also be determined from 
acoustic data acquired by the receiver. Setting up a seismic survey generally involves 
the step of acquiring acoustic data at the receivers, to verify that the receivers are 
located at their desired survey locations. This preliminary acoustic data can be used in 
the determination of the receiver heading by a method of fee invention. The receiver 
heading determined from the acoustic data can be used in the processing of seismic data 
subsequently acquired by the receiver. 

It is also possible to acquire seismic or acoustic data specifically for the purpose of 
determining the receiver heading. 
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As noted above, in a real receiver one of the horizontal sensing elements may measure 
the component of the wavefield more accurately than the other horizontal sensing 
elements so a calibration filter a(©)is generated to calibrate one horizontal sensing 
element against the other. Uncertainties in the calibration of the horizontal sensing 
elements introduce an ambiguity between the receiver heading and the required 

* 

calibration filter a(o>). In such a case, the cost function Ei should be minimised for data 
acquired over a number of shots covering a specific range of azimuths, rather than for 
data acquired at a single shot. Care should however be taken to minimise problems that 
might arise owing to azimuthally varying calibration filters, by selecting data for shots 
covering only a limited range of azimuths at a time. 

Once the sensor heading has been determined it is then possible to calibrate the 
measured pressure against the horizontal components of the particle velocity. This is 
done in a similar way to the determination of the sensor heading but, instead of having 
the sensor heading as an unknown variable, frequency dependent filters a(<o) and b(G>) 
are convolved with the x and y components of the particle velocity. Since the sensor 
heading is known it is possible to solve for a(a>) and b(<o), and this allows the x and y 
components of the particle velocity to be calibrated against the pressure. 



Any information about the calibration filters for the horizontal components obtained 
here is, in principle, available for quality control purposes during the acquisition stage. 
It can also be used during the subsequent processing of the acquired seismic data. 



The spatial derivatives in horizontal directions of the pressure can be obtained explicitly 
either by comparing data acquired at one receiver for two different source positions 
(source side derivatives), or by comparing data acquired at adjacent receivers for a 
single shot (receiver side derivatives), provided the data have been acquired with 
sufficiently small source or receiver spacing, respectively. 

Figures 1 and 2 shows possible seismic surveying arrangements for carrying out the 
method of the invention. Figure 1 shows a seismic surveying arrangement having a 
plurality of seismic sources (three sources 1,1,1 are shown in Figure 1), in this example 



WO 02/073239 



PCT/IB02/01897 



13 

towed by a vessel 2. The resulting pattern of shot points is not linear, as would be the 
case for a single source, but is a two-dimensional array, and this allows the horizontal 
derivatives of the pressure to be estimated on the source side so that the receiver 4 used 
in the seismic surveying arrangement of Figure 1 may be a conventional 4-component 
receiver. 



The seismic surveying arrangement of Figure 2 has only a single seismic source so that 
the horizontal derivatives of the pressure must be estimated on the receiver side, for 
example by using a receiver of the type described below with reference to Figure 4. 

The receiver 4 in Figures 1 and 2 may be disposed on the sea-bed, or it could be towed 
through the water. 

The method of the present invention requires only the phase differences between 
different recordings (that is, between recording of the same event at different spatially 
separated receivers or between recordings at one receiver of similar events from 
different, spatially separated sources). The phase differences can be obtained by 
determining the lag of the maximum of the cross-correlation of these recordings. The 
phase difference <pf thus obtained is equal to the product of the slowness p 3 - of the 
incoming wave and the source or receiver spacing Ai, where the subscript i denotes the 
direction (x or y) of the source or receiver separation. Assuming plane waves, the 
spatial gradients can be related to the temporal derivatives in the following way: 

Inserting the above in equation (S), and eliminating the time-derivatives allows us to 
write an alternative formulation of the cost function Ei, namely: 
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<p x v 



x y 



Ajc 



Ay J 



Up-sampling of the data prior to cross-correlation will increase the accuracy of the 
estimated phase-differences. 

In principle, in this embodiment the sensor heading could be determined by calculating 
the quantity /A x - #5, v yre c /Ay) and finding the orientation of the receiver axes 

that makes this quantity zero. In practice, as explained above, it is more convenient to 
find the orientation of the receiver axes from the square of this quantity. 

A conventional 4-component seismic receiver contains geophones for recording the 
three components of the elastic wavefield (these are generally embodied as a single 
multi-component geophone) and a hydrophone for recording the pressure. Since the 
pressure is recorded at only a single location it is not possible to obtain the spatial 
derivatives of the pressure on the receiver side if such a receiver is used. Such a 
conventional receiver could be used in the surveying arrangement of Figure 1 in which 
the derivatives of the pressure are determined on the source side. 

Once seismic data has been acquired in a seismic survey that uses a seismic surveying 
arrangement as shown in Figure 1 or 2, the receiver heading can be determined from the 
acquired seismic data using a method as described above. The determination of the 
receiver heading is conveniently done by processing a portion of the acquired seismic 
data, for example the time window containing the direct arrival. Once the receiver 
heading has been determined, some or all of the acquired seismic data can be processed 
in any conventional manner taking account of the determined receiver heading. 

Figure 4 is a schematic illustration of a seismic receiver that is suitable for obtaining the 
spatial derivatives of the pressure on the receiver side, and so is suitable for use in the 
surveying arrangement of Figure 2. The receiver 6 of Figure 4 has a multi-component 
sensor 7 for recording three orthogonal components of the particle displacement, 
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particle velocity or particle acceleration (in principle the multi-component sensor could 
record components of higher order time-derivatives of the particle displacement), and 
also has at least three hydrophones (or other suitable pressure sensors). Three 
hydrophones 8,9,10 are shown in Figure 4, but in principle more than three 
hydrophones could be provided. The multi-component sensor 7 and the hydrophones 
8,9,10 are disposed in a housing 11. 

The method of determining the receiver heading from the acquired data requires 
knowledge of only the x- and y~components of the particle displacement (or of the 
particle velocity etc), so the sensor 7 could in principle measure only the x- and y- 
components of the particle displacement (or of the particle velocity etc). Preferably, 
however, the sensor 7 is a three-component sensor so that the vertical component of the 
particle displacement (or of the particle velocity etc) can be acquired. The multi- 
component sensor 7 is preferably a multi-component geophone although, in principle, a 
plurality of single-component geophones could be used instead of a multi-component 
geophone. 

Figure 4 also shows the support cable 5. The receiver 6 may be coupled to the support 
cable 5 as in conventional receiver cable systems, or it may be decoupled from the 
support cable. Electrical connections from the support cable 5 to the sensor 7 and the 
hydrophones 8,9,10 have been omitted from Figure 4 for clarity. 

Figure 4 shows only one receiver 6, but in practice a large number of receivers would be 
provided on the support cable 5, at intervals along the length of the cable. 

The hydrophones 8,9,10 of the seismic receiver 6 are distributed over the horizontal 
plane in a non-collinear manner, so that the hydrophone distribution has a finite extent 
in both the Xrec and y m dimensions. The hydrophones preferably lie in a common 
horizontal plane. This hydrophone distribution makes it possible to obtain the 
derivatives of the pressure in the x Tec and y^c directions, by comparing the pressure 
values recorded at two or more of the hydrophones. For example, an estimate of P 
can be obtained by subtracting the pressure recorded at the hydrophone 10 from the 



WO 02/073239 



PCT/IB02/01897 



16 

pressure recorded at the hydrophone 8, and dividing by the distance between the two 
hydrophones. 

The accuracy of the estimate of the derivatives of pressure improves as the separation 
between the hydrophones is reduced. The separation between neighbouring 
hydrophones is therefore preferably of the order of centimetres, in which case the 
receiver 6 can be said to have a spatially dense areal distribution of hydrophones. 

Figures 1 and 2 show the receiver 4 in a marine seismic survey. A receiver of the 
present invention is not limited to use in marine seismic surveying, however but could 
also be used, for example, in a land-based seismic survey or deployed in a borehole. In 
addition to use in the method of the present invention, a seismic receiver of the 
invention is also advantageous for other purposes in hydrocarbon exploration, including 
noise attenuation and providing quality control measures for seismic data. 

An example of the method of the invention will now be illustrated by means of 
synthetic seismic data generated with a reflectivity code according to the method 
described by B. L. N. Kennett, in "Seismic wave propagation in stratified media" 
Cambridge University Press, Cambridge, England (1983). These synthetic seismic data 
were simulated for a seismic surveying arrangement of the type shown in Figure 2 in 

I 

which the pressure derivatives are obtained on the receiver side. The data input to the 
method consists of a time-window containing the direct arrival only. The simulation 
was performed for a receiver oriented with its Xrec- and y^c-axes coincident with the 
global x- and y-axes (i.e., for a = 0). The synthetic seismic data produced in the 
simulation included the pressure, the spatial derivatives of the pressure along the x^- 
and y Tec -axes, and the components of the particle velocity along the x^ and y^-axes. 

The middle frame of Figure 3(a) shows a cross-section through the cost function Ei over 
an azimuth range of -45° < 9 < +45°. The darkest areas of Figure 3(a) represent areas 
where the cost function is lowest The cost function Ei is calculated for noise-free 
synthetic seismic data, and is calculated over the range -10° < a < +10°. Since the 
simulation was carried out on the assumption that the receiver axes were aligned with 
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the global axes, the minimum of the cost function Ei should occur at a - 0. The central 
frame of Figure 3(a) shows that the minim um of the cost function does indeed occur at 
a = 0, over the entire azimuth range studied. The sensor heading is well defined and 
can be estimated accurately for all azimuths of the incoming wave. 

The upper and lower frames of Figure 3(a) illustrate the effect of incorrect calibration of 
the sensing elements for recording the horizontal components of the particle velocity. 
The simulations of these frames were performed using the same noise-free synthetic 
data as used to generate the central frame. The effect of incorrectly calibrated sensing 
elements were simulated by applying a calibration filter to the in-line horizontal 
component of the measured particle velocity (Vx). A simple calibration filter was used 
in the simulation, consisting of multiplying the in-line horizontal component of the 
measured particle velocity by a scalar factor sex before computing E* . The top and 
bottom frames of Figure 3(a) show cross-sections through Ei obtained from synthetic 
data in which v x has been multiplied by a scalar factor sc^ = 0.8 (upper frame) and by a 
scalar factor sc x =1.2 (lower frame). It will be seen that the minimum of Ei does not 
occur at zero (except at 0 = 0) in the upper and lower frames - incorrect calibration of 
the sensing elements leads to an error with a cosine-type behaviour in the estimated 
sensor heading. 

The estimated orientation of the sensor should be identical for all azimuths of the 
incoming wave. Thus, if the calculation of the cost function Ei should lead to results of 
the form shown in the upper or lower frame of Figure 3(a) and so show that the 
calibration of the sensing elements is incorrect, it is possible to use this fact to 
determine the calibration factor required to calibrate the sensing elements correctly. 
This may be done by finding the calibration factor that minimises the variation of the 
estimated sensor heading with azimuth, and this can be done by finding the calibration 
factor that minimises the variance, or some other parameter indicative of the variation, 
of the estimated sensor heading with azimuth. That is, the correct calibration factor is 
the calibration factor that minimises a second cost function E 2 : 

£ 2 = var(min(^ 1 )) . 
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Figure 3(b) shows the second cost function E2 as a function of sex, for the data of Figure 
3(a). As expected, Figure 3(b) shows that the required calibration factor in the 
simulation of Figure 3(a) is sc* = 1.0. Figure 3(b) shows that the desired calibration 
factor is well constrained and can be estimated accurately. 

In the example of Figures 3(a) and 3(b) the effect of incorrectly calibrated sensing 
elements has been simulated by a calibration filter that is simply a constant factor sc x . 
In a practical seismic survey the calibration filter is likely to be more complicated than 
this. However, the correct calibration filter can still be found by determining the 
variation of the estimated sensor heading with azimuth. The correct calibration filter is 
the calibration filter that minimises the variation of the estimated sensor heading with 
azimuth. 

Figures 3(c) and 3(e) correspond generally to Figure 3(a), but were obtained for 
synthetic seismic data in which the partial derivatives of the pressure and the 
components of the particle velocity contain white noise. The results of Figure 3(c) were 
obtained for synthetic seismic data having a signal-to-noise ratio of 20dB and Figure 
3(e) were obtained for synthetic seismic data having a signal-to-noise ratio of lOdB, 
Figures 3(d) and 3(f) show the second cost function E2 for the data of Figures 3(c) and 
3(e) respectively. In all these numerical examples, the sensor heading could be 
estimated well within an accuracy of 0.5 degrees. The form of the second cost function 
E2 is generally unaffected by the addition of noise to the synthetic data, although the 
minimum value of E2 occurs slightly away from scx= 1.0 in Figure 3(f)- 

The methods of calibrating the two horizontal sensing elements proposed by Gaiser 
(supra) and Bagaini et al (supra) are both based upon minimisation of energy transverse 
to the radial-vertical plane, where the radial direction is defined by the source and 
receiver co-ordinates. These methods require that the source and receiver co-ordinates 
are known, so that the radial direction can be determined, but this will not always be the 
case. 
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The present invention can be applied to these prior art calibration methods. The present 
invention allows the direction of the azimuth (6) of the seismic wave angle to be 
determined, and this angle defines the radial direction. Thus the invention allows the 
radial direction for each shot to be determined separately, thus avoiding the possible 
inaccuracies associated with the use of estimated source and receiver co-ordinates. 

The present invention has been described above with reference to marine seismic 
surveying methods in which the receivers are deployed on a sea-bed cable. The present 
invention is not, however, limited to marine seismic surveying. The method of 
determining the orientation of a sensor can in principle be applied to any multi- 
component seismic receiver such as, for example, seismic receiver deployed on land, a 
seismic receiver that is towed through water, and a seismic receiver deployed in a 
borehole. The method of calibrating the horizontal sensing elements of a seismic 
receiver can, in principle, also be applied to any seismic receiver. Moreover, the 
seismic receiver of the invention is not limited to use in OBC (ocean-bottom cable) 
seismic surveying arrangements but can be used in any conventional seismic surveying 
arrangement. 

In the example of the invention described above, the sensor heading is determined using 
spatial derivatives of the pressure in two different horizontal directions and using two 
horizontal components of the particle velocity* However, the invention is not limited to 
using components of the particle velocity. In its broadest form, the invention 
determines the sensor heading from information about at least one of the pressure 

* 

measured at the receiver and the particle displacement, or a derivative thereof measured 
at the receiver. In the above example, the first cost function Ei contains components of 
the particle velocity (that is, components of the time-derivative of the particle 
- displacement). However, the components of the particle velocity in the first cost 
function Ej could be replaced by components of the particle displacement, to produce an 
alternative first cost function E'i 
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where u xrec , are the components of particle displacement along the receiver's x- and 
y-axes. The sensor heading may be determined by finding the orientation of the 
receiver x- and y-axes that minimises the cost function E'i. 

Alternatively, the components of the particle velocity in the first cost function Ei could 
be replaced by components of the particle acceleration (that is, the second time- 
derivative of the particle displacement), to produce a further alternative first cost 
function E-'\ 

where a xreCi are the components of particle acceleration along the receiver's x- and 
y-axes. The sensor heading may be determined by finding the orientation of the 
receiver x- and y-axes that minimises the cost function E n i. 

In principle, the components of the particle velocity in the first cost function Ei could be 
replaced by components of the third or higher order time-derivative of the particle 
displacement. 

* 

Where the sensor heading is determined using the particle displacement or acceleration 
(or a higher time-derivative of the particle displacement), it is possible to calibrate the 
measured pressure against the components of the particle displacement, acceleration or 
higher time-derivative* This can be done using a calibration method similar to the 
method described above of calibrating the measured pressure against the horizontal 
components of the particle velocity once the sensor heading has been determined. That 
is, frequency dependent filters a(co) and b(co) are convolved with the x and y 
components of the particle displacement or acceleration (or higher time-derivative of 
the particle displacement) and, since the sensor heading is known, it is possible to solve 
for a(<x>) and b(co) and so allow the x and y components of the particle displacement or 
acceleration (or of a higher time-derivative of the particle displacement) to be calibrated 
against the pressure. 
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In the preferred embodiments described above, the receiver heading has been 
deteiraiaed using horizontal spatial derivatives of the pressure at the receiver. In 
principle, horizontal spatial derivatives of the particle displacement at the receiver could 
be used in place of the spatial derivatives of the pressure. If this is done, the sensor 
heading may be found using cost functions analogous to the cost functions Ei, Ei 1 and 
Ei" above, but with the terms d yrec P and dxreeP replaced by d yre <M and d xrBC u. In principle 
the spatial derivatives of the particle displacement may be determined on the source 
side, for example by using a two-dimensional array of sources as in Figure 1, or on the 
source side using a receiver analogous to that shown in Figure 4 but having a two- 
dimensional array of particle displacement sensors. 

Figure 5 is a schematic block diagram of an apparatus 12 for performing a method 
according to the present invention. The apparatus comprises comprises a programmable 
data processor 13 with a program memory 14, for instance in the form of a read only 
memory ROM, storing a program for controlling the data processor 13 to perform, for 
example, a method as described hereinabove. The system further comprises non- 
volatile read/write memory 15 for storing, for example, any data which must be retained 
in the absence of power supply. A "working" or "scratchpad" memory for the data 
processor is provided by a random access memory (RAM) 1 6. An input interface 17 is 
provided, for instance for receiving commands and data. An output interface 1 8 is 
provided, for instance for displaying information relating to the progress and result of 
the method. Data for processing may be supplied via the input interface 17 or may 
optionally be provided by a machine-readable store 19. 

The program for operating the system and for performing a method described 
hereinbefore is stored in the program memory 14, which may be embodied as a semi- 
conductor memory, for instance of the well-known ROM type. However, the program 
may be stored in any other suitable storage medium, such as magnetic data carrier 14a 
(such as a "floppy disc") or CD-ROM 14b. 
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CLAIMS 



1 . A method of detemiining the orientation of a seismic receiver, the method 
comprising: acquiring seismic or acoustic data at a receiver, the data including at least 
one of pressure and the particle displacement or a time-derivative of the particle 
displacement; and determining the orientation of the receiver from at least one of the 
pressure measured at the receiver and the particle displacement or time-derivative of the 
particle displacement measured at the receiver. 



2. A method as claimed in claim 1 wherein the step of determining the orientation 
of the receiver comprises determining the orientation of the receiver in a horizontal 
plane. 



3. A method as claimed in claim 1 or 2 wherein the step of determining the 
receiver orientation comprises: determining the spatial derivatives of the pressure at the 
receiver in first and second horizontal directions, the first direction being different from 
the second direction* 



4. A method as claimed in claim 1 or 2 wherein the step of determining the 
receiver orientation comprises: determining the spatial derivatives of the particle 
displacement at the receiver or a time derivative thereof at the receiver in first and 
second horizontal directions, the first direction being different from the second 

■ 

direction. 



5. A method as claimed in claim 3 or 4 wherein the step of determining the 
receiver orientation comprises: determining components of the particle displacement at 
the receiver in the first and second horizontal directions. 



6. A method as claimed in claim 3 or 4 wherein the step of determining the 
receiver orientation comprises: determining components of the particle velocity at the 
receiver in the first and second horizontal directions. 
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7. A method as claimed in claim 3 or 4 wherein the step of determining the 
receiver orientation comprises: detemiiriing components of the particle acceleration at 
the receiver in the first and second horizontal directions. 

8. A method as claimed in claim 6 when dependent from claim 3 wherein the step 
of determining the receiver orientation comprises detennining the quantity v^ec dy rec P - 
Vyrec ^xrecPy where Vxrec and v yrec are the components along the receiver's x-axis and y- 
axis respectively of the particle velocity measured at the receiver, and d^cP and d^P 
are the derivatives along the receiver's x-axis and y-axis respectively of the pressure 
measured at the receiver. 

9. A method as claimed in claim 8 and further comprising the step of detennining 
an orientation of the receiver's x- and y-axes that give a value of approximately zero for 

10. A method as claimed in claim 8 or 9 and comprising the step of determining an 
orientation of the receiver's x- and y-axes that minimises the function 

11. A method as claimed in claim 1 or 2 wherein the step of detennining fee 
receiver orientation comprises determining the quantity v^ec /A x - ft v ync /A>.) where 
v xrec and v^cc are the components along the receiver's x-axis and y-axis respectively of 
the particle velocity measured at the receiver, is a phase difference, and A* is the 
source or receiver spacing in the direction. 

12. A method as claimed in any preceding claim and further comprising the step of 
detennining the variation of the receiver orientation with azimuth. 



13 . A method as claimed in claim 12 and comprising the step of determining a 
calibration filter than minimises the variation of the receiver orientation with azimuth. 
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14. A seismic receiver comprising three or more pressure sensing elements, the 
positions of the pressure sensing elements not being collinear. 

15. A seismic receiver as claimed in claim 14 wherein the separation between any 
two of the pressure sensing elements is of the order of centimetres. 

16. A seismic receiver as claimed in claim 14 or 15 and comprising sensing 
elements for measuring two horizontal components of particle velocity. 

17. A seismic receiver as claimed in claim 16 and further comprising a sensing 
element for measuring the vertical component of particle velocity. 

18. A seismic receiver as claimed in claim 16 or 17 wherein the sensing elements 
for measuring the components of particle velocity comprise a multi-component 
geophone. 

19. A method of seismic surveying comprising: emitting seismic or acoustic energy 
at a first location; and acquiring seismic or acoustic data at a receiver disposed at a 
second location different from the first location; wherein the receiver is a receiver as 
defined in any of claims 14 to 18. 

20. A method of seismic surveying as claimed in claim 19 wherein the method is a 
method of marine seismic surveying and the receiver is disposed on the sea-bed. 

21. A method of seismic surveying as claimed in claim 19 or 20 and comprising the 
further step of processing a first portion of the acquired data according to a method as 
defined in any of claim 1 to 13 thereby to determine the orientation of the receiver. 

22. A method of seismic surveying as claimed in claim 21 and comprising the 
further step of processing a second portion of the acquired data taking account of the 
determined orientation of the receiver. 
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23. An apparatus for determining the orientation of a seismic receiver from seismic 
or acoustic data acquired at the receiver, the data including at least one of pressure and 
the particle displacement or a time-derivative of the particle displacement; the apparatus 
comprising means for determining the orientation of the receiver from at least one of the 
pressure measured at the receiver and the particle displacement or time-derivative of the 
particle displacement measured at the receiver. 

24. An apparatus as claimed in claim 23 and comprising means for determining the 
orientation of the receiver in a horizontal plane. 

25. An apparatus as claimed in claim 23 or 24 and comprising means for 
determining the spatial derivatives of the pressure at the receiver in first and second 
horizontal directions, the first direction being different from the second direction. 

26. An apparatus as claimed in claim 23 or 24 and comprising means for 
determining the spatial derivatives of the particle displacement at the receiver or a time 
derivative thereof at the receiver in first and second horizontal directions, the first 
direction being different from the second direction. 

27. An apparatus as claimed in claim 25 or 26 and comprising means for 
determining components of the particle displacement at the receiver in the first and 
second horizontal directions. 

28. An apparatus as claimed in claim 25 or 26 and comprising means for 
determining components of the particle velocity at the receiver in the first and second 
horizontal directions. 

29. An apparatus as claimed in claim 25 or 26 and comprising means for 
determining components of the particle acceleration at the receiver in the first and 
second horizontal directions. 
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30. An apparatus as claimed in claim 28 and comprising means for determining the 
quantity v xrec dyrecP - v yrec d xr ecP<> where v xr€C and v yrec are the components along the 
receiver's x-axis and y-axis respectively of the particle velocity measured at the 
receiver, and d xr€c P and dyreJP are the derivatives along the receiver's x-axis and y-axis 
respectively of the pressure measured at the receiver. 

31. An apparatus as claimed in claim 30 and comprising means for determining an 
orientation of the receiver's x- and y-axes that give a value of approximately zero for 



32. An apparatus as claimed in claim 30 or 31 and comprising means for 
determining an orientation of the receiver's x- and y-axes that minimises the function 

#1 = (Vxree dyncF — Vyrec dxrecP) 2 - 

* 

33. An apparatus as claimed in claim 23 or 24 and comprising means for 
determining the quantity v xrec /Ax - v yrec /Ay) where v mc and y^ c are the 
components along the receiver's x-axis and y-axis respectively of the particle velocity 
measured at the receiver, ^ is a phase difference, and At is the source or receiver 
spacing in the z ,th direction. 

34. An apparatus as claimed in any of claims 23 to 33 and comprising a 
programmable data processor. 

< 

35. A storage medium containing a program for a data processor of an apparatus as 
claimed in claim 34. 
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